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Abstract. We study the distribution of low-redshift Lya 
clouds in a CDM model using numerical simulations in- 
^vQ eluding photoionization and cooling of the baryonic com- 
ponent. The ionizing background is found to be efficient 
^ enough to keep most of the gas warm (T = 1-5 10 4 K) 
and to prevent most of it from collapsing. In this scena- 
t-H rio, the numerous low redshift (z < 0.5) neutral hydrogen 
Lya lines recently revealed by HST observations of QSOs 



originate in the warm gas that traces the potential wells 



of filamentary structures defined by the dark matter. The 
"■^ predicted line number density is consistent with observa- 
*^ tions. 

II About 25% of the Lya lines are found to arise in exten- 
O ded envelopes of luminous galaxies or groups of galaxies 
H (M > 3 10 11 M ) of radius of the order of 1.5 Mpc. The 



t/3 



^ exact number, however, is model dependent. The rema- 
inder arises in gas independent of any galaxy and at a 
distance of 2 to 7 Mpc from the first neighbour galaxy. 



association of the Lya lines with galaxies. They find 11 
galaxies coincident in redshift with an absorption system 
of which nine are Lya-only systems and two have asso- 
ciated CIV absorptions. They report that (i) most of the 
z ~ 0.3 luminous galaxies are surrounded by gaseous en- 
velopes of radius ~160/ij~ 1 kpc and (ii) the fraction of 
absorption clouds lying at a distance smaller than about 
150/ij~oo kpc from a luminous galaxy is at least 0.35±0.10 
and likely 0.65±0.18. However most of the Lya systems 
with associated galaxy have W T > 0.5 A , which transla- 
tes into Y(HI) > 5 10 16 cm" 2 for b = 20 km s" 1 (see their 
Figs. 23 and 26), whereas the bulk of the Lya lines has 
W T < 0.3 A (Morris et al. 1991). Indeed very recent HST 
observations of bright nearby AGNs by Stocke et al. (1995) 
have revealed one weak (W T = 36 mA) Lya line with no 
galaxy within 6 Mpc. The question of whether weak ab- 
sorptions are associated with galaxies or are completely 
unrelated to them is still open (Morris et al. 1993). 
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1. Introduction 

A large number of Lya clouds have been detected very re- 
cently at low redshift by HST (Morris et al. 1991, Bahcall 
et al. 1991, 1993). The observed number density of lines 
significantly exceeds the simple low-redshift extrapolation 
of the high-redshift z-dependence. This observation could 
be explained by a nonzero cosmological constant. Howe- 
ver, a significant fraction of the low-redshift population 
could represent relics of a population of clouds which, at 
high redshift, were associated with galaxies. Indeed, Lan- 
zetta et al. (1994) have determined the redshift of field 
galaxies around QSOs observed by HST to search for any 
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Galaxies are strongly clustered at low z. If the ab- 
sorbing gas is associated with galaxies, the absorption 
lines are expected to cluster on scales smaller than 
1000 km s _1 . The ID correlation should have an excess on 
those scales. The issue is unclear however. On the theo- 
retical side, since the size of the Lya structures may be 
larger than 10% of the mean separation between absor- 
bers, the spatial correlation function is diluted (Bajtlik 
1995). The observational facts are not clearly established 
as well. At high redshift (z > 2), early investigations have 
concluded that Lya lines do not cluster on scales lar- 
ger than 50kms _1 . However, new high resolution data 
(Cristiani 1995) show that the clustering properties de- 
pend on the column density threshold of the sample and 
a signal is detected in the ID correlation function of the 
gas with Y(HI) > 10 14 cm -2 for velocities smaller than 
200 km s _1 . This together with the very recent detection 
of weak CIV lines in Lya systems with column densities 
N(m) ~ 10 14 5 cm" 2 (Tytler 1995), argues for part of the 
Lya forest being associated with galaxies. Given the above 
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difficulties, the number of Lya lines detected at low reds- 
hift by HST is not large enough to allow a reliable statisti- 
cal study of clustering at low redshift. However, detection 
of clumps of Lya lines spread over about 1000 km s _1 has 
been reported by Boksenberg (1995). 

Since observational evidence for direct association of 
most Lya absorptions with bright galaxies is not unequi- 
vocal, we suggest here that they arise in gas associated 
more generally with the potential wells of rich filamentary 
structures that are ubiquitous in gravitational evolution 
dominated by dark matter. To obtain a consistent picture 
within the framework of models for large-scale structure, 
we study the spatial distribution of Lya gas and its time 
evolution using N-body simulations. In this paper we fo- 
cus our attention on the number density of Lya lines and 
the spatial distribution of the Lya gas at low redshift, in 
particular on its possible association with very extended 
disks or halos of galaxies. The properties of these systems 
at high redshift and their time evolution will be addressed 
in a forthcoming paper. 

2. Simulations 

We have improved and adapted the particle-mesh (PM) 
code described in detail by Kates et al. (1991; hereafter 
KKK) and Klypin & Kates (1991; hereafter KK) to in- 
clude the effects of photoionization. Let us briefly review 
the procedure, emphasizing those features relevant to the 
present investigation: 

The PM code simulates the kinematics of the dark 
matter and the thermodynamics of baryons. Here, we work 
in a (25 Mpc) 3 box (unless otherwise stated, h = .5 is as- 
sumed) with 128 3 particles on a 256 3 periodic grid. The 
positions, velocities, "temperatures", and gas densities as- 
sociated with particles are followed starting from a COBE- 
normalized CDM spectrum and utilizing the Zel'dovich 
approximation until it breaks down. 

The baryonic material is approximated as following 
trajectories of the dark matter with constant £lj = .1. 
This approximation is quite accurate in all regions prior 
to the first pancakes (or "shocks"). After shocking, the ac- 
curacy of the constant-Slj approximation is of course de- 
graded in regions where secondary shocking and/or very 
efficient cooling take place (e.g., in clusters and in the cen- 
tres of galaxies, see KK and KKK), but accuracy in these 
regions is not critical here. The conditions for maintai- 
ning constant S!j are well satisfied in Lya clouds, at least 
within the mini-halo picture (Rees 1986): The potential 
wells of mini-halos must be deep enough to prevent the 
gas from escaping, but still shallow enough to keep the 
central matter density below levels permitting strong coo- 
ling and subsequent star formation to occur. 

As perturbations grow, the first objects start to col- 
lapse, producing shocks in the gas. At a shock, the tem- 
perature assigned to the gas in a particle with velocity v 
entering a region with local velocity U is kT = / uh"^h(v — 



U) 2 /3, where fin is the molecular weight and mjj the mass 
of the hydrogen atom. We take primordial abundances 
throughout: Although as discussed above at least part of 
the Lya forest may contain metals, abundances cannot 
be so high as to significantly influence the quantities of 
interest here. 

After assignment of temperature, we follow the ther- 
mal history by integrating the energy equation along par- 
ticle trajectories, which takes the form 



dT < n 



T dn H hm 1 



dt 



(Atot — r p hot) 



(1) 



where -y=5/3 and «h is the number density of hydro- 
gen atoms. Here, A to t is the total cooling rate (see KK 
and KKK) including radiative and dielectronic recombi- 
nation, bremstrahlung, collisional excitation (Black 1981) 
and Compton cooling (the latter is the most important 
process at z > 5). 

The heating rate due to photoionization by the ionizing 
background (i.e., from QSOs and galaxies) is denoted in 
Eq.(l) by Tphot^ The heating rate is computed as the sum 
of the contributions by HI, Hel and Hell with 



r H i = F e- Tm Gen m 



where 
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for neutral hydrogen and similar relations for helium. 
Here, a is the cross-section for photo-ionization, and 
e is the mean kinetic energy per photo-ionization (va- 
lues are taken from Black 1981). The ionizing flux F v 
is modelled by a power law of index -1. F a is the ioni- 
zing flux at 13.6 eV and depends on redshift. We take 



F = /4tt10- 



Hz" 



2 at z ~ 2.5, where / 



is a dimensionless parameter with values varying in the 
range 0.2-2. Bajtlik et al. (1988) argue that / may be 
greater than one. On the other hand, from detailed analy- 
sis of metal line systems, Petitjean et al. (1992) found that 
this may well be an upper limit depending on the shape of 
the ionizing spectrum. A relatively high value / = 0.9 is 
required here in order to reproduce the redshift evolution 
of the line number density (Miicket et al. 1995). 

At redshift z > 5, essentially nothing is known about 
the physical state of the intergalactic medium and the 
ionizing UV background. In fact, as discussed in KKK 
and KK, due to the efficiency of Compton cooling at 
high redshift, the CDM model (and presumably all models 
with substantial small-scale power) requires some rehea- 
ting mechanism. Otherwise, most particles shocked prior 
to z = 5 would cool rapidly, implying excessive star forma- 
tion. Here we consider photo-ionization as such a mecha- 
nism. The redshift dependence of the ionizing background 
has been studied by several authors (Bechtold et al. 1987, 
Miralda-Escude & Ostriker 1990) and has been shown to 



be very sensitive to what are the dominant sources of pho- 
tons. It is clear however that the collapsed regions are the 
places where the photoionizing photons originate. 

To treat this question consistently, we take the ionizing 
flux intensity proportional to the rate at which material 
cools below 5000 K in the simulation. This simulates an io- 
nizing flux, produced as a consequence of gas collapsing in 
the dense regions, and then playing the role of a regulator 
preventing most of the gas in intermediate density regions 
from collapsing. The redshift dependence is thus consi- 
stently given by the simulation; the only free parameter 
here is the normalization: We assume / = 1 at z = 5. The 
resulting redshift dependence is found to be a power law 
(l+z) 13 with j3 = -2 and 2 for z > 5 and z < 2.5 respec- 
tively. Between z = 2.5 and 5, the flux is approximately 
constant. This is very close to what is expected in Models 
6 and 8 of Miralda-Escude and Ostriker (1990). 

In Eq. (2), r is a characteristic optical depth in the 
cloud. To estimate it, we model each particle as an in- 
dependent cloud of constant density taken as the local 
density given by the code. The radius is thus known from 
the mass of the particle, and r is the optical depth in the 
centre. This model would tend to overestimate the num- 
ber of strong lines, since gravity makes the cloud more 
centrally concentrated; the number obtained here should 
thus be regarded as an upper limit (see Section 3). The 
density ratio n(HI)/n(H) is computed self-consistently, ta- 
king into account photoionization and collisional ioniza- 
tion. Although thermal equilibrium is certainly not re- 
ached for i!h < 10 -4 cm -3 (see eq. 1), ionization equi- 
librium is a good approximation even for low densities 
(Duncan et al. 1991). 

When material cools below 5000 K, small, dense mole- 
cular clouds may form. Therefore the corresponding par- 
ticules are assigned the label "cooled" which is kept for 
the remainder of the simulation. In order to take into ac- 
count density variations through the cloud and subsequent 
evolution, a detailed hydrodynamic description would be 
required (see e.g. Cen et al. 1994). However specific stu- 
dies have shown that this is most important for the high- 
density regions of the clouds and that the low-density re- 
gions relevant to our study are only marginally affected 
(Murray & Petitjean in preparation). 



3. Results 

Fig. 1 shows the projected spatial distribution at z = 
of the dark matter particles in a 25x25 Mpc 2 slice of 
depth 2 Mpc. The structures seen form the usual net- 
work of filaments connected by nodes where galaxies and 
groups of galaxies are to be found. The neutral hydro- 
gen column density is computed through the slice, and 
contours of constant neutral hydrogen column density 
A(HI) = 10 13 cm -2 are shown. It can be seen that the 
Lya gas traces the gravitational potential wells. 
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Fig. 1. Projected space distribution of cold dark mat- 
ter particles in a slice of size 25x25x2 Mpc 3 at redshift 
z = 0. Contours for constant neutral hydrogen column density 
JV(HI) = 10 13 cm" 2 are overplotted. 



At high redshifts (2 < z < 3), the number density 
dN/dz of Lya systems is a strongly decreasing func- 
tion of time, dN/dz oc (l+z) 7 . Using high spectral reso- 
lution data, Rauch et al. (1992) found j = 2.1±0.5 for 
logA(HI) > 13.75 and in the redshift range 2.0 < z < 3.5. 
^From a compilation of low resolution data, Lu et al. 
(1991) derived 7 = 2.75±0.29 for lines with w T > 0.36 A. 
Extrapolating these observations would imply a very low 
number density at low redshift. Recent observations by 
HST have shown that the actual observed number of li- 
nes is five to ten times larger that what such a simple 
extrapolation would suggest. The number density of lines 
per unit redshift is 100±25 for A(HI) > 10 13 cm" 2 (Mor- 
ris et al. 1991) and 15±2 for W > 0.32 A (Bahcall et al. 
1993) which corresponds roughly (depending also on 6) to 
A(HI) > 10 14 cm -2 . In our model the number density of 
lines with A(HI) > 10 13 cm" 2 is 105 at z = 0. The ability 
to reproduce this number, for a value / = 1 within ob- 
servational constraints on the ionizing flux, represents an 
important success of the scenario. The number density of 
lines with A(HI) > 10 14 cm -2 is about 20 in the model. 

The idea that bright galaxies are surrounded by exten- 
ded envelopes is clearly verified in our simulations. Warm 
gas is present in large potential wells up to several hun- 
dreds of kiloparsecs from the centre (see Fig. 1). 

Observational evidence for the association of Lya gas 
with halos of galaxies was discussed in the Introduction. 
In order to quantify the "association" predicted by CDM, 
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Fig. 2. Histogram of the distance from Lya clouds to the near- 
est halo of mass M > 3 10 11 M 

we first identify all dark-matter halos in the (25 Mpc) 3 si- 
mulation box with M > 3 10 11 Mq, obtaining 13 of them. 
This mass range and spatial density lend credence to the 
assumption that the halos identified indeed represent indi- 
vidual "bright galaxies" (or groups of galaxies). We then 
compute the distance r of any cell within which Lya gas 
is found with 7V(HI) > 10 14 cm -2 to the nearest "bright 
galaxy" (the number of cells is proportional to the vo- 
lume occupied by the warm gas). A histogram of r values 
is shown in Fig. 2. The distribution is essentially bimo- 
dal. First, there is a peak at r < 1.5 Mpc corresponding 
to an excess of Lya gas around the galaxies themselves. 
This implies, as is apparent in Fig. 1, that bright gala- 
xies are indeed surrounded by halos of radius of the order 
of 1 Mpc. Secondly there is a broad but well-defined re- 
gime out to about 7 Mpc where the distribution drops off 
sharply which is consistent with a predicted characteristic 
mean filament separation of about 14 Mpc (Doroshkevitch 
& Turchaninov, in preparation). In this regime, the Lya 
gas is not related to any "bright galaxy". The volume oc- 
cupied by the gas lying at a distance greater than 1 Mpc 
to any bright galaxy is twice as large as the volume of 
the galaxy halos. If we assume, as a first approximation, 
that the halos around bright galaxies are spherical and 
the distribution along filaments is cylindrical, the ratio of 
the projected surfaces of both structures, S c /S s , is pro- 
portional to the ratio of their volumes multiplied by 
where R s is the radius of the sphere and r c the radius of 
the cylinder. Taking R s ~ 1.5 Mpc and r c ~ 0.8 Mpc gi- 
ves S c /S s ~ 3 and leads to the conclusion that about 25 % 
of Lya absorptions should arise within 1.5 Mpc of a bright 
galaxy. This result is of course sensitive to the operational 
definition of a "bright galaxy", a long-standing problem 
in N-body simulations. It is also sensitive to the local des- 
cription of the Lya gas. Both questions will be addressed 
in a forthcoming paper. 



Testing the picture presented here - that the Lya fo- 
rest traces the filamentary network of the dark matter - 
will require new observations. In particular high-quality 
data should be obtained for the Lya forest in QSO pairs 
with very different separations to probe scales from a few 
tens of kiloparsecs to several megaparsecs. The probability 
distribution of coincidences at different scales would pro- 
vide strong constraints on the model. On the other hand, 
a detailed investigation of the metal content of the Lya 
forest at high redshift using 10m class telescopes and at 
low redshift using shift and stack methods are thus clearly 
needed. 
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